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ABSTRACT 

Aims. Recent ALMA observations identified one of the most massive star-forming cores yet observed in the Milky Way; SDC335- 
MMl, within the infrared dark cloud SDC335.579-0.292. Along with an accompanying core MM2, SDC335 appears to be in the early 
stages of its star formation process. In this paper we aim to constrain the properties of the stars forming within these two massive 
millimetre sources. 

Methods. Observations of SDC335 at 6, 8, 23 and 25 GHz were made with the Australia Telescope Compact Array. We report the 
results of these continuum measurements, which combined with archival data, allow us to build and analyse the spectral energy 
distributions (SEDs) of the compact sources in SDC335. 

Results. Three Hyper Compact Hii regions within SDC335 are identified, two of which are within the MMl core. For each HCHii 
region, a free-free emission curve is fit to the data allowing the derivation of the sources’ emission measure, ionising photon flux and 
electron density. Using these physical properties we assign each HCHii region a zero-age main sequence spectral type, finding two 
protostars with characteristics of spectral type B1.5 and one with a lower limit of B1-B1.5. 

Ancillary data from infrared to mm wavelength are used to construct free-free component subtracted SEDs for the mm-cores, allowing 
calculation of the bolometric luminosities and revision of the previous gas mass estimates. 

Conclusions. The measured luminosities for the two mm-cores are lower than expected from accreting sources displaying character¬ 
istics of the ZAMS spectral type assigned to them. The protostars are still actively accreting, suggesting that a mechanism is limiting 
the accretion luminosity, we present the case for two different mechanisms capable of causing this. Finally, using the ZAMS mass 
values as lower limit constraints, a final stellar population for SDC335 was synthesised finding SDC335 is likely to be in the process 
of forming a stellar cluster comparable to the Trapezium Cluster and NGC6334 I(N). 
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1. Introduction 


The massive infrared d ark cloud (IRDC) SDC335.579-0.292 
( [Peretto & Fuller||2009| ) has been noted as an object of interest 
due to its large mass but low levels of radiative feedback for the 
protostars inhabiting it. The lack of the disruptive effects of ra- 
diative feedback allow SDC335 to be used to observe the initial 
. . conditions of star formation. 

^ SDC335 shows an interesting filamentary structure (Figure 
filamentary arms radiating outward from a denser 
central region. Within this central region reside two infrared 
^ bright cores, MMl and MM2. The cores are separated by 20.2" 
equivalent to 0.32pc at the distance of SDC335 w hich through- 


Peretto et al. 


out th is work i s taken as 3.25(^[j -^^)kpc (following 
|2013| using the |Reid et al.|2009| model). 

Bot h cores are associate d with 6.7 GHz class II methanol 
masers ( [Caswell et al.|201l] ), with two methanol masers associ¬ 
ated with MMl and a single maser associated with MM2. The 
class II methanol maser is a unique tracer of massive star for - 
mation ( jMinier et al.||2Q03| jXu et al.||2008t [Breen et ar||2Q13| ), 
meaning that SDC335 is harbouring massive protostellar ob- 
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jects. However, prior to this work no free-free emission at 6 cm 
had b een detected to a sensitivity limit of 0.2mJy ( [Garay et~ar 
2002 [ ), indicating that these sources are indeed young and in a 
pre-ionising phase of their evolution. 

SDC335 wa s the target of ALM A Cycle 0 observations with 
ALMA Band 3 ( [Peretto et al.[2013[ ). ALMA observed the 3.2mm 
dust continuum as well as spectral line transitions for CH 3 OH 
(13-12) and N 2 H^ (1-0). Using these data it was shown that 
MMl has a mass of ~500 Mq making it one of the most mas¬ 
sive star forming cores observed in the Galaxy. MM2 was found 
to have a mass of ~50 M(^ The ALMA data supplemented by 
Spitzer, Herschel and Mopra data also showed that the whole 
of SDC335 was undergoing global collapse, the converging net¬ 
work of filaments feeding the central region of the cloud with 
pristine gas. 

In this paper we present follow up observations of SDC335 
made using the Australia Telescope Compact Array in the radio 
( 6 , 8 , 23 and 25 GHz) and use these new data to better constrain 


^ We here state approximate masses for MMl and MM2 as a calcu¬ 
lation error has been found in [Peretto et al.| ( [2013| ) since publication 
leading to ~16% smaller masses. This issue will be addressed in a forth¬ 
coming corrigendum. 
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Fig. 1. Spitzer three colour image of the IRDC SDC335. Colour scale: 
red, green and blue GLIMPSE images at 8.0, 4.5 and 3.6//m respec¬ 
tively. White contours: ALMA 3.2mm continuum, with contours at 3, 
10, 25, 50, 75 and 90% of the peak flux. Green ★, Class II CH 3 OH 
masers. The yellow dashed lines outline the filaments of SDC335, fol¬ 
lowing |Peretto ( |2013| ) . 


Table 1. Observing characteristics of interferometric studies of SDC335 
used in this study. 


Telescope 

Freq. 

[GHz] 

Continuum 
Sensitivity 
[mJy bm“^] 

Synthesised 
Beam 
["X "] 

ALMA 

93.7 

0.400 

5.6X4.0 

ATCA 

25.0 

0.023 

1.9X1.0 

ATCA 

23.0 

0.026 

2 .1x1.1 

ATCA 

8.0 

0.008 

2.3X1.0 

ATCA 

6.0 

0.013 

2.9X1.4 


the masses of the stars forming in MMl and MM2. In addition, 
we discuss the implications of the presence of three newly iden¬ 
tified HCHii regions associated with these millimetre sources. 


2, Observations 

2.1. ATCA 

The primary source of data for this work is the Australia Tele¬ 
scope Compact Array (ATCA), from two dedicated observing 
runs under project code C2835. 

Centimetre wavelength (6 and 8 GHz) observations were 
made using the ATCA over 13 hours on 9^^/10^^ June 2013. 
During these observations the ATCA was in the ‘6C’ array con¬ 
figuration, giving a wv-range from ~2 to 140kA. These observa¬ 
tions consisted purely of continuum observations with the CABB 
system ( [Wilson et al.||201l] ) set up to observe a 2 GHz band¬ 
width centred at both 6 and 8 GHz simultaneously. PKS1934- 
638 was used for flux and bandpass calibration and PKS1646-50 
for phase calibration. 

Millimetre wavelength (23 and 25 GHz) observations were 
taken in a single 11 hour session on 13^^ September 2013, with 
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Fig. 2. PACS 160//m image of SDC335, overlaid with ALMA 3.2mm 
data (contours). It is clear that MMl and MM2 are blended in this image 
with emission from MMl possibly affecting photometry for MM2. The 
ALMA synthesised beam shape is seen in the lower left inset. 


the ATCA in the H.5A’ array configuration. During processing, 
data from the longest ATCA baseline (all antenna 6 baselines) 
was flagged out to allow better matching of t/v-range to the 6 and 
8 GHz data; the millimetre wavelength data therefore have a uv- 
range from ~9 to 127kT. The CABB system was configured to 
observe two 2 GHz continuum bands centred at 23 and 25 GHz 
plus seven ‘zoom’ bands to cover six ammonia transitions and 
the water maser at 22 GHz. Only the continuum are discussed in 
this paper, the ammonia and water maser data will appear in a 
later work. The flux and phase calibrators were the same as the 6 
and 8 GHz observations with PKS 1253-055 used for bandpass 
calibration. 


These ATCA data were reduced using the data reduction 
package MIRIAD with standard ATNF reduction strategies. The 
task MFCLEA10was used to CLEAN the image owing to the large 
bandwidths of the data. We refer the reader to the MIRIAD user 
guid^for more information. 


The beam and sensitivity characteristics of these observa¬ 
tions are given in Table including the jPeretto et al.j p013| ) 
ALMA observations of this object for comparison. Integrated 
flux density measurements from these observations were ob¬ 
tained by fltting Gaussians to each source using the MIRIAD 
IMF IT task, with the exception of MM 1^2 at 6 GHz where the 
slight blending of the two MMl sources led to unreliable gaus- 
sian fits, as such values for this source were measured “by hand” 
using a defined polygonal area in CASA ( jMcMullin et al.|2007] ) 
(leading to the increased uncertainty in the 6 GHz integrated flux 
density value). The integrated flux density values are shown in 
Tabled 


^ http ://www. atnf.csiro.au/computing/software/miriad/doc/mfclean.html 
^ http ://www. atnf.csiro.au/computing/software/miriad/userguide/ 
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Table 2. Observed results of the SDC335 milimetre and protostellar cores used in this study. 


Instrument 

Ereq. 

[Hz] 

A 

[m] 

Resolution^ 

["] 


^ int 

[Jy] 


MMl 

MMla 

MMl/? 

MM2 

PACS 

4.28 (12) 

70.0 (-6) 

5.0 

792 ± 158 

- 

- 

604 + 121 

PACS 

1.87 (12) 

160.0 (-6) 

12.0 

1456 ± 582.0 

- 

- 

961 + 384 

LABOCA 

3.45(11) 

870.0 (-6) 

18.0 

8.92 ± 5.34 

- 

- 

1.45 + 0.87 

ALMA 

9.37 (10) 

3.20 (-3) 

5.6 X 4.0 

1.01 ±0.10 (-1) 

- 

- 

1.20 + 0.20 (-2) 

ATCA 

4.59 (10) 

6.54 (-3) 

3.2 X 0.8 

6.85 + 1.76 (-3) 


- 

<0.92 (-4)* 

ATCA 

4.32 (10) 

6.95 (-3) 

3.5 X 0.9 

6.56 + 1.37 (-3) 

3.14 + 0.91 (-3) 

- 

<0.89 (-4)* 

ATCA 

2.50 (10) 

1.20 (-2) 

1.9 X 1.0 

2.21 + 0.32 (-3) 

1.41 +0.17 (-3) 

7.12+ 0.12 (-4) 

7.17 + 0.99 (-4) 

ATCA 

2.30 (10) 

1.30 (-2) 

2.1 X 1.1 

2.48 + 0.35 (-3) 

1.76 + 0.22 (-3) 

6.02+ 0.12 (-4) 

7.24+ 1.10 (-4) 

ATCA 

8.00 (09) 

3.75 (-2) 

2.3 X 1.0 

6.07 + 0.98 (-4) 

3.02 + 0.39 (-4) 

3.00 + 0.36 (-4) 

3.59 + 0.42 (-4) 

ATCA 

6.00 (09) 

5.00 (-2) 

2.9 X 1.4 

3.33 + 0.52 (-4) 

1.60 + 0.36 (-4) 

2.25 + 0.25 (-4) 

2.72 + 0.34 (-4) 


t, for interferometers resolution given is the synthesised beam major and minor axes. 
★, 3cr upper limits on the none detection of MM2. 

Columns 2, 3 and 5-8 values in brackets are the order of magnitude in base 10. 


2.2. Additional data 

In addition to the ATCA observations listed above, data for 
SDC335 was available at 7mm from ATCA (Avison, Cunning¬ 
ham & Fuller, in prep. Note only MMl is d etected in these data) 
and ALMA at 3.2mm (|Peretto et al.||2013|). Archival data from 
Spitze^MlPS and GLIMPSE), MSX, and APEX (LABOCA) 
have also been used within this work. 

We also made use of Herschel data from the "Herschel in¬ 
frared Galactic Plane Survey" (Hi-GAL, Molinari et al. 2010), 
which has mapped the whole Galactic Plane covering the wave¬ 
length rang e 70 < T < 5 00 pm using both t he PACS ( [Poglitsch] 
|et ^.||2010| ) and SPIRE ( [Griffin et al.||2010| ) photometry instru- 
ment s in parallel mode on board of the Herschel space observa¬ 
tory ( jPilbratt et al.|2010| . We limited our analysis to the PACS 
70 and 160 pm data, which have the highest spatial resolution 

5” and ^ 12” respectively). 

The source extraction and photometry of the Herschel data 
has been carried out with Hyper, an enhanced aperture photom¬ 
etry code specifically designed to account for high background 
variability and source crowding of the Galactic Plane data ( jTraf-j 
jicante et al.|2014| ). The background is evaluated locally for each 
source assuming different polynomial orders (from zero to the 
fifth order) to model the background fluctuations. The polyno¬ 
mial background-subtracted map with the lowest rms is taken as 
reference to evaluate the source flux, and the equivalent polyno¬ 
mial background order as the best fit to describe the background 
fluctuations in the Hi-GAL maps. 

In the case of overlapping sources, as is the case of MMl and 
MM2, prior to evaluating each source flux in the aperture region 
the algorithm models all the sources with a multi-Gaussian fit 
and subtracts the model of the companions. Within SDC335 the 
proximity of the brighter MMl core to MM2 at the resolution of 
Herschel means that the two cores are blended in both the 1.9 
and 4.3THz (160 and IQpm) Herschel maps (see Figure]^. It is 
the data points from these maps which most constrain the SED fit 
(e.g. Figure discussed later in this work. The flux estimation 
for each core can be affected by the flux of the overlapping com¬ 
panion or by a coarse background estimation. Although Hyper 
can estimate fluxes with a precision of few percent even in com- 


^ Data from Spitzer is shown within Figures and for illustrative 
purposes. 


plex background regions and in the presence of heavily blended 
sources ( jTraficante et al.||2014| ), we decided to be conservative 
and assume a flux uncertainty of +20% of the photometric Her¬ 
schel flux values. 

The source fluxes extracted with Hyper have been corrected 
for aperture corrections. The integrated fluxes from this extrac¬ 
tion, Sint, are shown in Table 

3. Results and analysis 


SDC335 as viewed in the radio/mm of our observations is seen 
in Figure Comparing our ATCA observations with the ALMA 
data we observe sources at the position of the two mm-cores, 
MMl and MM2, in SDC335. However, given the increased res¬ 
olution of our ATCA data we find that the MMl core com¬ 
prises two sources separated by ~3". We name these MMl a 
and MMl/? , with MM la being the more southerly of the pair 
(see Figure [^b). The two fragments are situated 2.7" from one 
another, a physical separation of -9000 AU at the distance of 
SDC335. MM2, the less massive mm-core in SDC335 remains 
a single source at this increased resolution. The position of each 
radio core is given in Table 

The red/blue contours and arrows in Figure 3 represent the 
5% of peak emission boundary and direction of the HNC molec¬ 
ular outflow discovered in the ALMA Cycle 0 data (to be dis¬ 
cussed in full in Pineda et al. in prep), with the colours denot¬ 
ing red and blue shifted emission respectively. Note there is a 
slight kink in the red shifted outflow, thus the two arrows for this 
feature. Of interest to the work in this paper is the fact that the 
position of maximal overlap between the two outflow arms co¬ 
incides, to within 0.5", with the position of peak continuum flux 
from MM la. As such we assume that MM la is the source which 
drives the molecular outflow. 

In Figure we also show the positions of three maser 
species: Class I (44 GHz) and Class II (6.7 GHz) CH 3 OH 
(with positions from Avison, Cunningham & Fuller, in prep, and 
[Caswell et al.|(|2011|), re spectively) and 22 GHz H 2 O (positions 
from [Breen et al.[ ( 2010[ ) and the 22 GHz spectral data from our 
ATCA observations). Class II CH 3 OH and H 2 O masers are typ¬ 
ically associated with the central regions of protostellar objects. 
Class II CH 3 OH masers are seen to trace discs or disc-like fea- 
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O) _ 


(b) 


- -48°43'48.0" 


MMla 



MMlb 


O.OSpc J 


59.04s 58.80s 58.56a6h30m58.32s 

RA (J2000) 



- -48°43'36.0" 


57.60s 57.36s 57.12a.6h30m56.88s 

RA (J2000) 


Fig. 3, (a) The SDC335 IRDC as viewed in the radio/mm continuum. Colour-scale the 8 GHz (ATCA) overlaid with ATCA 25 GHz (magenta 
contours at 15, 20, 30, 50,75, 80 and 90% of the peak flux) and the ALMA Cycle 0 3.2mm continuum (grey contours at 2.5,10, 25, 50, 75 and 90% 
of the peak flux). The location of maser emission in SDC335 are indicated by x, • and -h for Class II (6.7 GHz) and Class I (44 GHz) CH 3 OH and 
H 2 O (22 GHz) masers, respectively. The red/blue contours and red/blue arrows represent the 5% of peak HNC emission boundary and direction 
of the HNC molecular outflow (Pineda et al., in prep) see main text for more detail). The synthesised beam shape for both the ALMA data (grey 
hatched) and ATCA 25 GHz (solid magenta) are shown in the lower left, (b) & (c) Zoom in of the central regions of MMl and MM2 respectively, 
showing 8 and 25 GHz continuum as per main image. 


tures when observed on milli-arcsecond scales (e.g. Bartkiewicz 
et al.|2009] ) and are uniquely associated with regions of massive 
star formation. H 2 O masers are typically associated with outflow 


regions in both high and low mass star forming regions (Fish 


2007), particularly the inner regions (-6200AU) thereof (Breen 
et al. |20I0[ and references therein) when the H 2 O and class II 


CH 3 OH masers are co-spatial. As such, each pairing of these 
maser species (one for each of the three radio continuum source 
in SDC335), is indicative of the position of a forming massive 
star. 

Class I CH 3 OH masers have been found to tra ce the larger 
scale molecular outflows of protostellar systems ( [Kurtz et~ar| 
2004[ l^ganowski et al.|2009| ). Two of the four class I CH 3 OH 
maser spots in SDC335 inhabit the blue outflow arm originat¬ 
ing from MMla. A third spot is adjacent to the HNC blue lobe 
observed near MM2. 

Given the assumed outflow origin, the angular offset between 
the outflow direction and a line connecting MMla and MMlb, 
and the presence of both a H 2 O and class IICH 3 OH maser asso¬ 
ciated with the peak continuum flux of MMla and b, we consider 
these two sources independent objects. 

Our analysis focuses on constraining the properties of 
MMla, MMl/? and MM2. 
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Table 3. The positions of the radio cores observed in SDC335. 


Source 

J2000 



Right Ascension 

Declination 


[ hh:mm:ss ] 


MMla 

16:30:58.765 

-48:43:54.00 

MMl/? 

16:30:58.638 

-48:43:51.70 

MM2 

16:30:57.291 

-48:43:40.21 


3.1. Free-free component 

The three compact Hii regions detected in our ATCA obser¬ 
vations are observed with sufficient resolution to create partial 
SEDs from 6 to 25 GHz, these can be seen in Figure (b) and 
(c) and Figure We calculate for each sources spectral index, 
a (where Sint oc v^), with Sint being the integrated flux of each 
source and v frequency, the a values are given in Table One 
can consider a as an indicator of the free-free emission turnover 
frequency. The turnover frequency of free-free emission being 
the frequency at which a source transitions from optically thick 
to optically thin, i.e. where the source optical depth equals 1. Be- 
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Table 4. Speetral index properties of the cores within SDC335. We do 
not include the a value between 23 and 25 GHz for sources MM la 
and MMIZ?, as it is possible the value is contaminated by some thermal 
emission (see Figure]^ (a)). 


y 

interval 

MMla 

a 

MMIZ? 

a 

MM2 

a 

6 - 8 GHz 

1.90+0.40 

1 . 01 + 0.12 

0.96±0.11 

8 - 23 GHz 

1.67+0.20 

0.66+0.13 

0.67±0.10 

23 - 25 GHz 

- 

- 

- 0 . 11 + 0.02 


low this frequency a free-free emission region is optically thick 
and a = 2, toward the turnover frequency a reduces and for op¬ 
tically thin sources a approaches -0.1 ( |Kurtz|2005] ). 

Each source is optically thick between 6 and 8 GHz. MM la 
appears to remain at least partially optically thick between 8 and 
23 GHz and the lower a values between these frequencies for 
MMIZ? and MM2 (compared to MM la) are also suggestive of 
partial optical thickness. However, these lower a values do allow 
for the turnover frequency to occur between these 8 and 23 GHz 
for both MMIZ? and MM2. 

MM2 is the only source for which there is insignificant con¬ 
tamination from thermal dust emission at 25 GHz, (contrast Fig¬ 
ure]^ and Figure [^, the a value between 23 and 25 GHz shows 
this Hii region has become optically thin at these frequencies. 

Based upon these spectral indices we consider the nature of 
the free-free emission, specifically considering the ambiguity be¬ 
tween compact Hii regions and collimated ionized jets at these 
resolutions. The steep spectral index of MM la bet ween 6-8 and 
8-23 GHz is typical of that of a compact Hii region ( |Kurtz|20Q5] ) 
and steeper than the expected value for collimated jet emission 
(Qf=0.25-l.l, [Reynold^ ( |1986| )). Thus we classify MMla as a hy¬ 
per compact Hii (HCHii) region. 

The spectral indices of MMIZ? and MM2 of ~1.0 and -0.6 
at 6-8 and 8-23 GHz respectively in both sources, fall within the 
chara cteristic range of both collimated jet emission ( [Reynolds] 
|1986| ) and photo-ionized compact Hii regions in the transition 
from optically thick to thin. In the case of MMIZ?, given its mis¬ 
alignment to the molecular outflow originating from MMla and 
centrally located (potentially disc tracing) CH 3 OH maser we as¬ 
sume MMIZ? is not part of a radio jet within the MMl mm-core 
and instead favour the interpretation that this source is also com¬ 
pact Hii region ionized by a separate source to MMla. We also 
assume the compact Hii region interpretation for MM2, based on 
the lack of a firm detection of a nearby strong collimated outfiow. 
However, within the scope of these data either the radio jet or 
compact Hii region interpretation may be valid for both MMIZ? 
and MM2. The ultimate result of this ambiguity is that the stel¬ 
lar mass limits calculated assuming a compact Hii region, ergo a 
photoionization origin for the radio emission, will not hold true 
for radio jets. However the association of MMIZ? and MM2 with 
class II CH 3 OH masers confirms that they are forming massive 
stars. 


3.1.1. Free-free fitting technique 

To further this analysis we fit the free-free component using the 
6 , 8 , 23 and 25 GHz ATCA data, implementing a non-linear least 
square regression algorithm. The fitted parameters are the free- 
free turnover frequency, v^, (where the optical depth, r = 1) and 
the physical diameter of the source, /. Setting r = 1 and fixing 
Te to the fiducial 10,000K (e.g |Spitzer|[r978t jSanchez-Mongej 


jet al.| 20 lT| ), emission measure, {EM), values for a given turnover 
frequency were calculated using the approximate formula for the 
optical depth of free-free emission ( jMezger & Henderson|1967 
and references therein): 


Tff ^ 0.08235 


GHz/ 


EM \ 
pc cm“^ / 


( 1 ) 


Based on this EM value the optical depth, Tff, following Equa¬ 
tion!^ was used in Equation 


5 ff.fit = nffB,{Te){l - e-^ff), (2) 

to fit the flux density of the data. Here ^// is the solid angle 
subtended by the source and By{Te) the Planck distribution value 
at a given frequency at 5//-//? was then fit to the observed 
fluxes with the best fit giving values for v^, EM and /. Note that 
the EM value calculated at the turnover frequency is independent 
of any knowledge of the sources physical dimensions meaning 
there is no degeneracy in the value of I which is used only in 
calculating the value of Elff. 


3.1.2. Fitted Properties 

The fitted properites for the three Hii regions seen in SDC335 
are given in Table Each object has fitted size of < 0.5mpc 
(^llOAU), a value much smaller than that of the upper size 
limit associated with the hyper-compact Hii region classification 
(<30mpc e.g |Kurtz|2005| ). Sources of similar size have bee n ob- 


served in high mass star forming regions, e.g. W3 IRS 5 (Wil- 
json et al]|2003| ) with sources dl/d2 and f both measuring de- 
convolved radii of 0.8mpc (^ 165AU) using the VLA at 1.3 and 


Q.7cm. [Zhu et al.| ( |2Q13| ) observed the HCHii region NGC7538 


IRSl with the SMA and CARMA (at 1.3 and 3.4mm) and re 
solved it into a compact component <270 AU (^1.3mpc) and an 
extended tail of 2000AU (^ 9mpc). Also, |Sanchez-Monge et al. 
(2011) in their VLA search for new HCHii regions (in a methanol 
maser selected sample) found six new HCHii candidates with de¬ 
rived sizes of between 1 and lOmpc (^206-2060 AU). 

The fitted optically thick/thin turn over frequency for the 
three SDC335 Hii regions are in agreement with the values ex¬ 
pected from the spectral indices calculated between each of our 
observing frequencies. MMla is the densest of the 3 Hii regions 
in SDC335 remaining optically thick to frequencies >20 GHz. 
MM 1/7 and MM2 exhibit similar properties to one another with 
fitted turnovers of - 11 GHz. From the fitted values for each 
of our compact radio sources it seems likely that SDC335 har¬ 
bours three HCHii regions. The fit to MMla is the least robust 
as the data for this source allows a rise in fiux beyond 25 GHz. 
However, the fitting technique limits the turnover frequency to 
be within t he freq uency range of the data. This issue is addressed 
further in §3.2.2| 


3.1.3. Derived properties 

For each source we calculated an electron number density, 
and the number density of Lyman a photons being emitted. The 
latter allowing the association of the three HCHii regions with a 
zero-age main-sequence (ZAMS) stellar spectral type. Assuming 
a cylindrical homogeneous Hii region EM is related to electron 

number density by EM =jl n^dl ^ nil. Calculated values for 
He are found in Table and again exhibit values associated with 
HCHii (> lO^cm'^). 
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Fig. 4. Spectral energy distribution of the MMl core on ~5"scales. a) The black o give the fluxes at each frequency, the blue x are the data used 
for the blackbody fit, note that at 23 and 25 GHz the free-free component has be subtracted. The black-body fit with and without the free-free 
component subtraction are given by the solid green and red dotted lines respectively. The magenta and cyan dot-dash lines and hexagons give the 
radio fluxes and follows the free-free fits for MMla and MMIZ? , respectively. The free-free emission data for MMla and MMIZ? can be seen 
separated in b) and c). 



= 7.6 X 10^'’ a-\v,Te) 




^ int,cont 

jy 


D 

kpc 


. 2.0 

photons S ^ 


(3) 


( |Zhu et al.|2013| ), where S int,cont is the radio continuum flux den¬ 
sity at a given frequency and D the distance to SDC335 (i.e. 
3.25kpc). For each source we calculate Niym using the 23 GHz 
continuum integrated flux density value. It is important to note 
that as MM la is optically thick beyond 25 GHz, use of the 23 
(or indeed 25) GHz will lead to an underestimate in Niym for this 
source meaning our calculated value is a lower limit . As previ- 
ously, we set re=10,000K, giving a(v,Te) = 0.966 ( [Mezger & 
|Henderson|1967| ). 


Observed Freq. [Hz] 

Fig. 5. Spectral energy distribution of the MM2 core. Symbols as Figure 

Eb. 


The Lyman continuum flux assuming optically thin emission 
is given by: 


The Niym values were then compared with the Lyman Con- 


tinuum values for ZAMS stars given in 

Davies et al.|(|2011|) and 

spectral types from (Mottram et al.|201 ] 

[|). From these values the 


protostellar cores in SDC335 appear to have currently formed a 
single B 1.5 (M^ = 9.0 M©) type star in the MM2 core, while the 
MMl core has formed at least two stars one of B 1.5 (MMl/?) 
spectral type and a second of at least spectral type Bl-B 1.5 
= 10.0 Mq) for MMla. These spectral classifications are based 
on the properties of the protostars as they are observed now were 
they at the ZAMS stage and not still actively accreting. 
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Table 5. Properties of the cores within SDC335 from the SED fitting. 


Object 

EM 

[pc cm“^] 

He 

[cm“^] 

Ihw 

[AU] 

Lym 

[photon s“^] 

Tbol 

[Lo] 

Td 

[K] 

p 

^gas 

[Mq] 

MMl 

- 

- 

- 

- 

1 . 8 x 10 ^ 

28.7(+2.3) 

1.97(+0.15) 

763C;®®) 

MMlp 

- 

- 

- 

- 

7.6x10^ 

30.0(±1.7) 

1.75(±0.11) 


MMla 

1.9(+0.5)x 10® 

2.1(+0.3)x 10® 

94.28(+5.7) 

>2.0(:^ 5) X 10^5 

- 

- 

- 

- 

MMl Up 

2.0{±0.4)x\Qp 

2.7(±a2)xl0® 

93.96(±7.1) 


- 

- 

- 

- 

MMIZ? 

5.0(±1.0)x 10*^ 

1 . 0 (± 0 . 1 )x 10 ® 

104.4(+5.2) 

6 . 8 (+' ®) X lO^'* 

- 

- 

- 

- 

MM2 

4.5(±0.5) X 10^ 

8.9(±0.5) X 10® 

115.2(+4.7) 

8.2C‘;®) X 10^4 

9.9x10^ 

29.6(+7.2) 

2.26(+0.52) 

173C“) 


t: ^Lym values are lower limits, see text. 

The values shown in rows MMlp and MMlUp are those which include the 43 GHz data point in the free-free emission SED fitting, see ^3.2.2| for 

more details. 


3.2. Spectral Energy Distribution of SDC335 protostellar 
cores 

Combining the multi-frequency data available (from this work 
and ancillary data described in |2.2| ), and listed in Table we 
create spectral energy distributions (SEDs) for the dust emission 
from MMl and MM2. These were created by first subtracting 
the best fit free-free emission from the thermal emission at all 
frequencies used in the fit (25 GHz to 4.3 THz (lOjim) for each 
source) and then fitting the SED to these corrected data. 

Eor the purpose of SED fitting MMl is treated as 2 objects 
(MMla and MMIZ?) at frequencies between 6 and 25 GHz and 
as a single object above these frequencies where the resolution 
is 5" or greater. As the respective resolutions of the data used 
in this fitting vary as a function of frequency, we conservatively 
compensate for potential inclusion of excess fiux in lower reso¬ 
lution data. To achieve this the uncertainties of the data points 
(presented in Table with beams greater than 6 " have been 
scaled assuming a density profile of r~^ for these sources (see 
[Peretto et al.|2013| Table 1.). 

The modified blackbody fitted is of the form: 


Sbb-fit = ^cBy(Td)(l-e-^^^) 


( 4 ) 


where By(Td) is the Planck distribution value at a given fre¬ 
quency, y, and a dust temperature T^. r^b is the optical depth 

of dust and has the form (—j , with vq as the frequency at which 


the optical depth is unity. Einally, is the solid angle subtended 
by the dust dominated region, for the purposes of fitting we fix 
Qc to the area sub tended from the Hersc hel data using the Hyper 
photometry code ( Traficante et al.|20T4] ) which has radius 11.4". 


The free parameters in our fitting routine are, vq and p, 
and are given priors of 25.OK, 4.3 THz and 2.0 respectively. The 
resulting best fit parameters are then used in a direct integration 
of the SED between 100 kHz and 100 THz (3km to 3yum) to 
calculate the source bolometric luminosity. We also refine the 
mass of the mm-cores MMl and MM2 using the 3.2mm ALMA 
fiux densities (see Table [^, at deconvolved s ize scales of 0.054pc 
and 0.057pc (following [Peretto et al.|2013] ), applying our fitted 
Tj, Vo and p in the mass flux relation. 


^ D^S^,2 

K 3 . 2 B 3 . 2 iTd) 


( 5 ) 


and using the dust opacity law = 0.1 x ( 0 3 '^^ ) g 


assuming a dust to gas ratio of 0.01 (|Hildebrand| 1983[|Beckwith 
|etal.|1990| >. 


3.2.1. SED fitting results 


Table presents the parameters derived from our SED fitting. 
Overall the values fitted and derived values from our SEDs are 
typical of those of a massive protostellar core in early stages of 


Star formation (e.g 

Rathborne et al.||2008[ [Beuther et al.||2010[ 

[Rathborne et al.|20] 

lOj. 


3.2.2. Posteriori fitting of MM1a 

As noted in §3. 1.1 [ the free-free emission from MM 1^2 is the least 
well constrained with the turnover potentially occurring above 
the 25 GHz data point used in the inital fit. To attempt to over¬ 
come this we use the SED fitted results for MMl (T=28.7K and 
P=1.91) following the above method. We then adjust the 43 GHz 
data point and refit the free-free emission from MM 1^2 including 
this data point. The model fit of MMlZ? is kept as constant within 
these fits. The fitting of free-free emission data and then the full 
(free-free subtracted) SED proceeds iteratively until best fit val¬ 
ues for MMlf 2 in free-free (/, Vt and EM) and new best fit values 
for MMl as a single object (T and P) are found. These new val¬ 
ues can be found in Tablej^labelled subscript ‘p’ and represented 
in italics. 

The new fitted properties for MMl and MMl a are consistent 
with the initial fitted values for each source, which may suggest 
a robust fit to these data or highlight the detrimental effect of the 
larger error associated with the 43 GHz data point compared to 
those from the primary data (6-25 GHz ATCA) for this work. 
As such we consider only the initial fit values throughout the 
remainder of this work. 


4. Discussion 


4.1. The evolutionary status of the ionising sources in 
SDC335 


In |3.1.3| we calculate a Lyman photon flux for each HCHii and 
from this associate a ZAMS spectral type to the source ionising 
each HCHii. We now consider the evolutionary status of the three 
ionising sources in SDC335 with res pect to the calculated Lboi 
values for the MMl and MM2 cores ( |3.2.1[ ) and the protostellar 
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Table 6. Calculated Luminosity values for MMla, MMIZ? and MM2 
using the assumed ZAMS properties of the ionising souree each. 


Source 

L* 

[i^o] 

^acc 

[i^o] 

Ltot,ZAMS 

[Le] 

MMla 

5.48x10^ 

2 . 22 x 10 =’ 

2.27x105 

MMIZ? 

4.10x10^ 

2.16x105 

2.20x105 

MM2 

4.35x10^ 

2.17x105 

2.21x105 


evolutionary tracers associated with each source compared to the 
characteristics that the ZAMS spectral type would apply. 

A star arrives at the Zero Age Main Sequence once hydro¬ 
gen burning has commenced. Within SDC335 we expect accre¬ 
tion to be ongoing, as evidenced by the presence molecular out¬ 
flow from the MMl mm-core (Pineda et al. in prep.), and each 
ionising source exhibiting a 22 GHz water maser ( [Breen et al 


( 2010| ); Avison et al. in prep.), a maser species associated with 


SDC335 filamentary structure. The global collapse of SDC335 
also supports the ongoing accretion of the protostars ionising the 
three detected HCHii. 

Accretion rates of between > and up to few x 


10 ^ M© yr ^ (e.g. 

Fuller et al. 

2005 

Churchwell et al. 2010| 

Klaassen et al.[2012 

IRygletal. 2013 ] 

Duarte-Cabral et al.|2013|) 


have been calculated for massive star precursors (e.g. HMPOs, 
HCHii), a greater value than that associated with their lower 

( Shu| [ 


mass counterparts 


5 X 10“'’ Mo yr- 


1977 


[McKee & Tan|[200^ ). Given such observed values we make 
the assu mption that the infa ll rate of SDC335 (2.5 + 1.0)xl0“^ 
Moyr“^ ( Peretto et al.|2013 ) continues unimpeded onto the pro¬ 
tostars. During the accretion phase of a young stellar object the 
total output luminosity is comprised of an accretion component 
and the intrinsic luminosity, Ltot = + Lacc with the accretion 

component taking the form Lacc = • If we assume the pro¬ 

tostars have the ZAMS radii, and luminosity, interpolated 
from Davies et al. (2011 Table 1) to our calculated NLym values 


and use the ZAMS mass-luminosity relation of oc we can 
calculate a mass for each protostar. With these values we then 
calculate Ltot for MMl and MM2 (in the case of MMl the value 
is the addition of calculated values for MM la and MMIZ?). In 
both cases Ltot exceeds the calculated by a factor ~ 20 , with 
Lacc the dominant luminosity (see Table [^. 

Initially we note that the ZAMS L^, without the additional 
luminosity contribution from Lacc, for both MMl and MM2 are 
within a factor ~ 2 of the measured Li^ou not a large discrep¬ 
ancy given the uncertainites in the data. However, the number 
of outflow indicators associated with the protostars (as noted 
above) excludes the scenario in which the SDC335 protostars 
have ceased accretion. Instead we consider two mechanisms by 
which the discrepancy between Ltot,zAMS and L^oi can be ad¬ 
dressed. 

First, we consider the models for massive protostellar evo- 
lution with high accre t ion ra tes ( > 1 O~^M 0 yr~^) presented by 
[Hosokawa & Qmukaij ( [2009[ ) and [Hosokawa et al.[ ( [201Q[ ) for 
spherical collapse and accretion via a disc respectively. Under 
both scenarios massive protostars transition through a number of 
stages before joining the ZAMS. Notably swelling at masses 6 
< < 10 M© and contracting again between 10 < < 30 M©, 

(for a constant M ~ 10“^M© yr“^). During the swelling phase the 


outflows from both low and high mass YSOs ( [Fish[[2007[). In 
addition to these, an Extended Green Object (EGO, CyganowsE[ 
et al.| 200 ^ , another outfiow tracer, is also present at the centre of 


protostellar radii swells rapidly from a few Rq to ~ 100 Rq whilst 
the intrinsic luminosity also increases. As a result of the swelling 
the effective temperature of the protostar is lower than that of 
a ZAMS star with the same luminosity, with temperatures low 
enough that there are insufficient photons with energies greater 
than 13.6eV to create a Hii region. 

Again following the M ~ 10“^ M© yr“^ case reported on 
for both accretion scenarios in [Hosokawa & Omukai[ ( [20()9| ) and 
Hosokawa et al.[ ( [2Q10[ ), after the swelling of the protostar a pe¬ 
riod of contraction commences whilst accretion continues. L^ 
becomes dominant component of Ltot over Lacc and R^ drops 
from ~100 7?© to of the order 10 7?© as the protostar increases 
in mass. It is during this phase that the protostellar effective tem¬ 
perature becomes high enough that the flux of ionising photons 
is sufficient to start the formation of a Hii region. 

Under this high accretion rate scenario and given that our 
protostars are displaying ionization but have Ljjoi < Ltot,zAMS 
we suggest the possibility that the protostars within SDC335 are 
in the contraction phase of their formation, putting their current 
masses each at > lOM© (see e.g. [Hosokawa & Qmukai[ ( [2009 


Figure 18)) with their eventual final ZAMS mass higher than 
this. It is interesting to note that our assumed mass accretion rate 
(2.5 + 1.0)xl0~^ Mn vr~^ approaches the limiting case of the 


Hosokawa et al. (2010) model at ~4xl0 M© yr ^ an accretion 


rate value at which a second stage of rapid radial expansion oc¬ 
curs halting contraction and steady accretion ceases. In this case 
the temperature again decreases owing to the radial expansion 
and no Hii regions are created. This provides a maximum limit 
on the accertion rates of these sources of < 4 x 10“^ M© 
( [Hosokawa et al.|2010[ ). 

The second mechanism which can address the luminosity 
discrepancy is that there is a decreased, decreasing or periodic 
accretion when compared to the used in our Lacc calc ulation. 
The value = 2.5 x 10“^M© yr“^ used was found by 


lo yr“ 


Peretto 


et al.[ p013[ ) and assumes spherical accretion onto the central 


mm-cores of SDC335 due to the global cloud collapse. It is pos¬ 
sible that the accretion rate onto the protostars is lower than into 
the cores, meaning that the cores are growing at a greater rate 
than the protostars thus Lacc is lower. Indeed, reducing the accre¬ 
tion rate onto the protostars by the factor 20 notes gives an M,. 
of 1.25 X 10“^M© yr“^ which is not unreasonable for a forming 
massive star. 

Decreasing or periodic accretion onto the protostars would 
lead to fluctuations in Lacc over time. These are the most diffi¬ 
cult scenarios to address in any quantitative manner as it requires 
both knowledge of the accretion on to the sources in a more di¬ 
rect way than our current data allow and preferably observation 
of this accretion or suitable proxies at more than one epoc h. Op¬ 
posing these accretion rate scenarios [Davies et al.[ ( [201 1[) found 
in their population synthesis based on the RMS data ( [Hoare et al. ' 
[2004[ ) that declining (or indeed constant) accretion rates as a 
function of time do not well fit the observed massive star popu¬ 
lation of the galaxy. There is currently no concensus on periodic 
accretion onto massive protostars within the literature. 

Our current data on SDC335 do not allow us to select the 
likely candidate from these scenarios as being the mechanism 
responsible for the lower than expected L^^/. Further observation 
of the SDC335 protostars targetted at constraining the accretion 
on smaller scales is required to resolve this issue. 

4.2. Cluster formation in SDC335 

Given that we have discovered three massive protostars within 
SDC335 it is possible to put limits on the eventual number of 
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stars the IRDC is capable of creating. We e mploy the power law 
form o f the IMF, N o c with a'=2.35 ( [Salpeter 1955 ) and 
a=23 ( |Kroupa|[2002| ) in the standard mass range of l-120Mo. 
With the criteria that 3 objects have masses above 9.7Mo, we 
find that SDC335 has the potential to form a minimum of ~54- 
60 stars (between 1 and I 2 OM 0 ). 

Continuing the |Kroupa| ( |2002| ) IMF down to a mass of 
O.O 8 M 0 SDC335 has the potential to form -1400 stars, with 
the vast majority of these (-96%) in the sub-solar mass regime. 
Using the average stellar mass from |Kroupa| ( |2002| ) Table 2 
this would give an approximate total mass in stars, of 

-285Mo, a small fraction (5%) of the total mass of the SDC335 
cloud ( 55 OO+ 8 OOM 0 ). Assuming a typical star formation ef- 
ficiency of -10 - 30% for embedded clusters ( |Lada & Lada 
|2003| ) one would expect a M^^tot from SDC335 in the range 
550-1650Mo. Keeping in mind the uncertainties on our mass es¬ 
timates, the discrepancy between this and our IMF calculated 
value could be the result of our underestimation of the true final 
masses of the three massive protostellar objects in SDC335. 


4.2.1. The future of SDC335 

For a young system such as SDC335, in the early stages of its 
star formation process it is interesting to make comparisons to 
systems with similar characteristics at different stages of evo¬ 
lution. One such system is the Trapezium Cluster (TC) at the 
centre of the Orion Nebula Cluster. This extensively studied rel¬ 
ative evolved system (see e.g. |Hillenbrand|1997t['Muench etaLl 
|2002[ ), has some obvious similarities with our you nger system 
Centered on O^C Ori, TC has a diameter of -0.3pc (Hillenbrand 


1997| ) comparable to the separation of the two SDC335 millime¬ 


tre cores (0.32pc). The eponymous Trapezium stars comprise 
fiye OB stars at the heart of the cluster similar to the three proto¬ 
stars exciting the HCHii regions described in this work. Extend¬ 
ing out to include the inner 3pc (and encompassing the TC) is 
the Orion Nebula cluster (ONC), a region of physical size com¬ 
parable to the whole of SDC335 (-2.4pc). 

Beyond these obyious similar ities, the ste l lar po pulation 
study of the ONC as reported by |Hillenbrand| \\991) detects 
-1600 optically yisible stars, with an expected additional 50% 
not detectedd due to their embedded nature. Of these sources the 
authors report on the 60% of stars for which photometry and 
spectral types are ayailable allowing these objects to be asso¬ 
ciated with a mass. In the mass range l-50Mo their sample is 
complete and contains 137 objects with 6 stars of mass greater 
than lOMo, (the Trapezium stars plus others at a larger distance 
from O^C Ori). We haye shown SDC335 has 3 stars of mass > 
IOMq and we expect 60 stars in a similar mass range to the 137 
stars in ONC, from this it would appear we can expect a similar 
population in SDC335 to that in the ONC to within a factor -2. 

It is also interesting to note that in SDC335 all three HCHii 
regions are found in the IRDCs central region (at the confluence 
of the filamentary arms). In the more eyolyed ONC Hillenbrand] 
( 1991) find eyidence for mass segregation with the most mas- 
siye stars preferentially located within projected radii of <0.3pc. 
In a following paper [Hillenbrand & Hartmann] ( |1998| ) find that 
this distribution of stars is unlikely to be the result of dynamical 
mass segregation since its timescale is larger than the age of the 
ONC stars. The authors instead fayour a primordial origin of the 
observed mass segregation. This interpretation is consistent with 
our SDC335 findings, the three most massiye sources lying at 
the bottom of the gas-dominated, grayitational potential well. 

A second comparable region is the NGC6334 I(N) proto¬ 
cluster which has been well studied at mm wayelengths, most 


recently by [Hunter et al.| ( |20T^ . NGC 6334 I(N) (hereafter sim¬ 
ply I(N)) is an IRDC at a distance of l.Skpc and displays seyen 
compact millimetre continuum sources within a projected diam- 
eter of 0.1 pc at 2"resolution. Further to this, in their study [Hunter 


[et al.[ ( [20T^ find 25 likely protostellar objects associated within 
a physical radius of 0.21pc at 1.3mm. 

There are four radio (5 GHz) continuum sources within the 
same 0.2Ipc radius area, three of which are labelled as proto¬ 
stellar candidates and two are directly associated with 1.3mm 
cores. The limited number of radio continuum detections in I(N) 
has an interesting corollary with SDC335, scaling the fiuxes of 
the I(N) radio continuum detections to the distance of SDC335 
we find that two of the four I(N) sources would be observed as 
5(T or greater in SDC335. Giyen this and the fact we do obserye 
three radio continuum sources in SDC335, if we assume I(N) and 
SDC335 haye similar protostellar population density profiles one 
would expect a similar number of compact mm sources to exist 
in SDC335. Howeyer, currently, the comparatiyely low resolu¬ 
tion of the [Peretto et al.[ ( |20T3] ) ALMA data at 3.2mm compared 
to that with of the 1.3mm obseryations of I(N) (coupled with the 
fact that these two mm frequencies will be tracing different re¬ 
gions around the protostellar populations therein) means we are 
limited in the our current knowledge of the mm-core population 
of SDC335 to make a quantatiye comparison. 

Qualitatiye comparisons of SDC335 to the NGC 6634 I(N) 
protocluster and the relatiyely eyolyed Trapezium cluster indi¬ 
cates that SDC335 is likely to form a similar population of stars. 
Further inyestigation of SDC335 with higher resolution at mm 
wayelength obseryations will allow the lower mass/non ionising 
core population to be detected making these comparisons more 
robust. 


5. Conclusions 


We haye presented the results of radio continuum obserya¬ 
tions conducted with the ATCA toward the infrared dark cloud 
SDC335. Within the two massiye protostellar cores observed by 
ALMA in SDC335 we haye obseryed three compact Hii regions. 
Fitting of the free-free emission of these sources indicates that 
they haye characteristics typical of Hyper Compact Hii regions. 
Based on the fitted and deriyed characteristics of these HCHii 
regions we find that the three protostars exciting them display 
characteristics of ZAMS stars of at minimum, spectral type B 1.5 
for two of the three and B 1-1.5 for the third. This can be con¬ 
sidered a lower limit of their characteristics as accretion is still 
ongoing within SDC335. 

With the use of ancillary data, free-free component sub¬ 
tracted SEDs were created allowing the calculation of the bolo- 
metric luminosities for the MMl (treated as a single object at 
frequencies aboye 25 GHz) and MM2 cores. We find yalues from 
our SED fits to each HCHii region and mm-core which are typi¬ 
cal of massiye protostellar sources. 


Our obseryed luminosities for the two cores are lower than 
would be expected from the ZAMS classification deriyed from 
the ionized flux. We interpret this result as a consequence of 
one of two mechanisms reducing the contribution of accretion 
to the total luminosity. First, the protostars may be at an early 
stage of accretion, still with an enlarged but contracting stellar 
radii after the ‘swollen core’ stage for high accretio n rate ob¬ 


jects (M=, > 10 ''^Mo yr ^ Hosokawa & Omukai 2009) indicating 


that the protostellar objects exciting our HCHii regions will ulti¬ 
mately become more massiye that their current ionising fiux in¬ 
dicates. Secondly the assumed accretion rate onto the protostars 
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within the mm-cores MMl and MM2 maybe lower than onto the 
cores themselves. 

Comparisons between SDC335 and both the realtively 
evolved Trapezium cluster and the evolving NGC 6334 I(N) 
protocluster have been drawn to highlight the similarities of 
SDC335 to these clusters and indicate the potential future of the 
cluster forming in SDC335. 

Use of future high resolution spectral line and dust contin¬ 
uum observations of SDC335 would allow detection of more 
precise knowledge of accretion rates of each object and detec¬ 
tion of the mass fragmentation and indeed segregation. These 
would provide more details of how SDC335 is forming a stellar 
cluster similar to the Trapezium Cluster and NGC 6334 I(N). 
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